/* Copyright (c) 2021, National University of Defense Technology. All rights reserved.
 *
 * Licensed under the Apache License, Version 2.0 (the "License");
 * you may not use this file except in compliance with the License.
 * You may obtain a copy of the License at
 * 
 *     http://www.apache.org/licenses/LICENSE-2.0
 * 
 * Unless required by applicable law or agreed to in writing, software
 * distributed under the License is distributed on an "AS IS" BASIS,
 * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
 * See the License for the specific language governing permissions and
 * limitations under the License.
 */

#include "coo_matvec.h"

namespace YHAMG
{

void COO_MatVec(const COO_Matrix& A, const Vector& x, const Vector& y)
{
	int n = A.size[0];
	int nnz = A.nnz;
	int* Ai = A.rowind;
	int* Aj = A.colind;
	double* Av = A.values;
	double* xv = x.values;
	double* yv = y.values;

	for (int i = 0; i < n; ++i)
		yv[i] = 0.0;

	for (int j = 0; j < nnz; ++j)
		yv[Ai[j]] += Av[j] * xv[Aj[j]]; 
}

}